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Abstract. We argue that isotropization and, consequently, thermalization of 
the system of gluons and quarks produced in an ultrarelativistic heavy ion 
collision does not follow from Feynman diagram analysis to any order in the 
coupling constant. We conclude that the apparent thermalization of quarks and 
gluons, leading to success of perfect fluid hydrodynamics in describing heavy 
ion collisions at RHIC, can only be attributed to the non-perturbative QCD 
effects not captured by Feynman diagrams. 

We proceed by modeling these non-pertrubative thermalization effects using 
viscous hydrodynamics. We point out that matching Color Glass Condensate 
inital conditions with viscous hydrodynamics leads to continuous evolution of 
all components of energy-momentum tensor and, unlike the case of ideal hy- 
drodynamics, does not give a discontinuity in the longitudinal pressure. An 
important consequence of such a matching is a relationship between the ther- 
malization time and shear viscosity: we observe that small viscosity leads to 
short thermalization time. 
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1. Introduction: Isotropization Versus Free Streaming 

The results presented here are mainly based on the work done in [HE], but also 
include new developments described in Sect. 4. 

Similar to the original Bjorken hydrodynamics approach [ 3 , let us consider a 
central high energy collision of two very large nuclei. For simplicity, here we will 
discuss the case where the distribution of particles is independent of space-time 
rapidity r] — (1/2) ln(a;+/j:_), where x± — {t ± z)/^/2. Since the nuclei are very 
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large the transverse coordinate dependence can also be neglected for most physical 
quantities, leaving only the dependence on the proper time t — y/2x^X- . For this 
geometry, one can show that the most general energy-momentum tensor can be 
written as (at z = 0) 




where the z-axis is taken along the beam direction, and x, y-axes are in the trans- 
verse direction. Applying the conservation of energy-momentum tensor condition 

df.T'"' = (2) 

to the energy-momentum tensor that gives Eq. at z = we obtain 

de _ e + P3 
dr T 

There are two interesting cases one can consider: 

(i) if P3 = longitudinal pressure vanishes and, due to Eq. we get 

1 



(3) 



r 



(4) 



such that the total energy E k er = const. This case is known as free 
streaming: the system expands freely without loosing any energy. 

(ii) if p3 = p the energy-momentum tensor in Eq. becomes isotropic. This is 
the case of ideal Bjorken hydrodynamics [ 3 . Eq. (|3J) with P3 = p was derived 
in [IS]. If combined with the ideal gas equation of state, e = 3p, it gives 



r4/3 



(5) 



or, for other equations of state, 

1 



with A > 0. (6) 

Eq. demonstrates that changes in the total energy E k, er (or, equivalently, 
deviations from e ~ 1/t scaling) are due to work done by the longitudinal pressure 
Pa. The classical initial conditions in the Color Glass Condensate (CGC) approach 
[0] yield the free streaming final state with p3 = 0. A thermalized quark-gluon 
plasma (QGP) is characterized by non-zero P3, leading to the energy density scaling 
as shown in Eq. ©. Therefore, below we will understand isotropization, which is 
the necessary condition for thermalization, as dynamical generation of non-zero 
longitudinal pressure ps ^ 0, or, equivalently, deviations from the scaling of Eq. 
leading to the scaling of Eq. © . 
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2. Formal Argument 

An extensive search of the diagrams which would bring in the desired deviations 
from the scaling of Eq. |0J carried out by the author did not yield any positive 
results: while many diagrams have contributions to e scaling as shown in Eq. 
such terms are always subleading additive corrections to the leading (at late times) 
terms scaling as shown in Eq. Q . In fact one can construct an argument demon- 
strating that the leading contribution to energy density from any-order diagrams 
scales as e ~ l/r. The argument is presented below. 

We begin by considering a gluon field generated by an arbitrary Feynman di- 
agram illustrated in Fig. ^ In d^A^^ — covariant gauge it can be written 
as 

where the function J^{k) denotes the rest of the diagram in Fig. ^ (the truncated 
part), which depends on the momenta of other outgoing gluons as well. Indeed 
gluon field can be defined as a simple function only in the classical case: the "field" 
in Eq. |71l should be thought of as a Feynman diagram in Fig. with one of the 
outgoing gluon lines being off mass-shell, i.e., a generalization of the classical field 
which we will need in calculating energy density [ ^ . (The expression in Eq. ((TJ 
can also be thought of as an operator equation.) 

Substituting Eq. Q into the expression for energy-momentum tensor 

Tf"" = l^-F'^f'P F'^'p + i {F^af^ , (8) 

averaging over the nuclear wave functions and employing the symmetries of the 
collision of two identical nuclei we obtain the energy density due to the gluon field 

m 



{2tt)» {k^ + ieka){k'^ + iek'f^) 
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where fi{k'^,k''^,k ■ k'jkx) is some unknown function (a "form- factor" ) and the 
elhpsis indicate addition of two more similar terms with different "form-factors" /2 
and /s. 

Rewriting each "form-factor" as 



k ■ fc', kr) = fi{k^ = 0, k'^ = 0, /c • fc' = 0, A:t)+ 
+ [/,(fc2, A: • k', kr) - hik" = 0, = 0, A: • fc' - 0, A:t)] 



(10) 



and using the fact that the square of truncated part of the diagram gives a cross 
section 

' ' r''{k)j^{-k) 



dPkdy 2(27r)3 



(11) 



we conclude that, keeping only the first term on the right hand side of Eq. (jlOfl for 
all "form-factors" in Eq. (O yields 



d^k 



dN 



d^k drj d^h 



krj^ 



[[Jl{kTT)f + [UkTT)f] « i j d^k 



dN 



d^k drj d^b 



where the last equality is valid for late proper times r. Since, as was shown in 
each factor of fc^, k'^ or fc • fc' gives a factor of 1/t, the terms in the square brackets 
of Eq. Hl()|l give a subleading (compared to Eq. p2(l ^ contribution to energy density 
at late times r and can be safely neglected.^ We have shown that any diagram 
and/or any set of diagrams contributing to gluon production cross section lead 
to energy density scaling as in Eq. i.e., that isotropization and, consequently, 
thermalization do not take place in perturbation theory analysis of the collisions. 

The main assumption of the argument presented above is the existence of mul- 
tiplicity of produced gluons dN/d^kdy, which is the essential assumption of QCD 
perturbation theory. This is what makes our argument perturbative. 



3. Physical Argument 

Now let us present a physical argument demonstrating the origin of the power of 
4/3 in Eq. (jSJ and explaining why it is impossible to achieve in perturbation theory 
[E]. Let us assume that thermalization does take place at some time Tth- If a 
gauge invariant time Tth exists, we can put the QCD coupling constant .9 = for all 
times r > Tth without violating gauge invariance. Bjorken hydrodynamics in the 
g = limit is governed by the ideal gas equation of state e — 3p, which leads to the 
energy density scaling as shown in Eq. (For small but non-zero g, Eq. |SJ) would 
get an o{g^) negative correction to 4/3: the expansion in g would still be around 

^We assume that fi's are continuous functions of their arguments. 
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the power of 4/3.) Due to Eq. the scaling of Eq. in the g = limit of 
Bjorkcn hydrodynamics means that ps and the gas of non-interacting particles 
is doing work in the longitudinal direction! What causes such a behavior of the 
system? The problem lies in the ideal gas equation of state, e = 3p, which assumes 
that the ideal gas is in contact with some external thermal bath. Such external 
thermal bath could be a background field or a box containing the gas: the ideal gas 
of non-interacting particles stays thermal through the interactions between the gas 
particles and the thermal bath. This is the only interaction allowed in the g — > 
limit and it is responsible for the work done by the non-interacting gas. Since there 
is no such external thermal bath in heavy ion collisions, the scaling of Eq. ((SJ is 
impossible to achieve at small coupling. 

Without the external thermal bath the particles in the gas would be just free 
streaming, giving the physically correct energy density scaling of Eq. Q). Of course, 
at a fixed time r hydro is not applicable in the .g = limit, since the mean free 
path of the particles would exceed the longitudinal size of the system. However, if 
thermalization does happen, for any fixed arbitrary small g, if we wait long enough 
hydrodynamics should become applicable, leading to the scaling arbitrary close 
to that of Eq. Q and doing work in the longitudinal direction which would be 
mostly due to contact with the non-existing external thermal bath. Therefore, we 
arrive at a contradiction, demonstrating that hydrodynamics is not achievable at 
small coupling. At large coupling, non-perturbative effects may mimic the external 
thermal bath, possibly leading to energy density scaling shown in Eq. J^Jl. 

4. Matching Color Glass Initial Conditions and Viscous Hy- 
drodynamics 

If thermalization does take place via some non-perturbative mechanism, under- 
standing its dynamics would be extremely hard. Instead, let us ask a pragmatic 
question: what is expected from a successful thermalization scenario? Apart from 
giving us correct dynamical mechanism for the generation of longitudinal pressure 
and the onset of isotropization, one may expect a thermalization scenario to yield us 
an estimate of the time when energy-momentum tensor would become symmetric, 
7)t!/ = diag{eiPjP,p} at z = (again we are considering a central collision of two 
heavy ions in the rapidity-independent approximation for produced particles). At 
this time ideal hydrodynamics would be initiated and it could be used to describe 
the subsequent evolution of the system. Unfortunately, as was argued in ideal 
hydrodynamics is probably unachievable in nature, since the shear viscosity r] never 
completely vanishes.^ Therefore, one has to use viscous hydrodynamics with the 
energy-momentum tensor given by 

Tf_,i, = {e + p) Uf, u„ -pgt^i, + ri ^V^u^ -|- V^u^ - ^ A^^ Vp u''^ , (13) 

^Following the standard convention we will use -q to denote shear viscosity in this chapter, 
which should not be confused with space-time rapidity -q used in the previous chapters. 
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where is the velocity profile, = (.g^,y — u^) and A^,y = g^i, — u^Ui,. In 
the boost-invariant case Eq. leads to 



One can see that viscosity corrections increase transverse pressure and decrease 
longitudinal pressure in the system. 

Now, if we follow [i5j and assume that ideal hydrodynamics is unachievable, we 
can try to address once again the question of what to expect from a thermalization 
scenario. Our earlier answer would not work anymore: thermalization can not give 
us a completely isotropic energy-momentum tensor, since such a tensor is impossible 
due to viscosity corrections. Therefore, thermalization dynamics can only bring the 
system to a somewhat isotropic state achieved at certain thermalization time, after 
which viscous hydrodynamics takes over. 

However, the advantage of viscous hydrodynamics energy-momentum tensor in 
Eq. I|14|l is that it can be continuously mapped onto the energy-momentum tensor 
of Color Glass initial conditions without including any additional thermalization 
dynamics. In the rapidity-independent case the energy-momentum tensor due to 
the classical fields in the Color Glass is given by [ 01 ^ 



We mark the energy density and pressure in Eq. H15|) with a prime to distinguish 
them from the appropriate components of the energy-momentum tensor in Eq. H14|l . 

Indeed matching the energy-momentum tensors in Eqs. (|15() and (|14|) is a 
rather crude approximation, which, by omitting thermalization dynamics, would 
only give us a lower bound on thermalization time. Alternatively, one can view 
viscous hydrodynamics (|14|l as a model for the non-perturbative thermalization 
effects: then, matching Eqs. (|15|) and H14|) would generate a non-perturbative 
thermalization scenario. 

Requiring that at some "thermalization" time tq the components of energy- 
momentum tensors in Eqs. (|15(l and H14(l are equal yields the following set of 
equations:'^ 



The energy density e' and transverse pressure p' of the CGC initial conditions 
determine, through Eqs. (|16|l . the initial values of energy density, pressure and 
viscosity of the QGP, along with the matching time tq. Indeed, the exact knowledge 
of QGP thermodynamics would, in principle, allow one to express e, p and rj as 

■^Continuous matching conditions in Eq. 1161 also lead to continuity of de/dr due to Eq. Isl . 




(14) 



T'^'' = diag{e'(r),p'(r),p'(r),0}. 



(15) 




(16) 
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functions of the system's temperature T. Then Eqs. would have only two 
unknowns - matching time tq and the corresponding temperature at the matching 
To, and would be over-constrained. However, in the case when Tj^ = on both 
sides of the matching, Eqs. (|16|1 have only two independent equations 

e(To) « e'(To) 

p(ro) - I ^ ^'^^ 

and are not over-constrained anymore. The tracelessness condition, = 0, is valid 
in CGC both at the classical level, and at the level of leading logarithmic small- a; 
evolution. On QGP side = is vaHd at very high temperatures, and also e = 3p 
appears to be a good approximation for the equation of state in the strong coupling 
regime. We conclude that solution of Eqs. H17(l in terms of tq and To is likely to 
satisfy the third equation in (fT?^ as well. 

It is impossible to solve Eqs. (|17ll without knowing the exact dependence of e 
and p on temperature. Instead, we will find an approximate solution by writing 

e{T) = 3piT) = ^nig)T\ (18) 

where n(g) is some function of QCD coupling constant g{T), which, at high tem- 
peratures, corresponding to small couplings, just counts the number of quark and 
gluon degrees of freedom. Following [ 5 we write for the shear viscosity 

71 ^ f{g'N,)N^T' (19) 

with f{g^ Nc) some function of the coupling. For the Color Glass energy density 
we write 

CfQ's 
as T 2tt'-^ 

where ce is some order 1 coefficient to be determined by explicit calculations and 
Qs is the saturation scale. Substituting Eqs. lfTH|) . ((T^ and if^ into Eqs. lfT7|l . and 
neglecting, perhaps unreasonably, the dependence of the coupling constant g{T) on 
the temperature T, we obtain the matching time 

^ [f{g^N,)N!]^/^ 120 ( ^\ ^ (21) 



>n{g) TT^ \ceCf 

which is a good estimate of thermalization time being the only relevant time scale. 

Let us analyze the thermalization time from Eq. H21(l . First of all, while, as we 
argued in the previous Sections, perturbative thermalization is dynamically impos- 
sible, we can still explore the small-coupling asymptotics of Eq. (|21|l to construct 
a lower bound on perturbative thermalization time. To that end we note that at 
small values of the coupling n{g) approaches a constant, while f{g^Nc) ^ 
Using this in Eq. H21|) . and dropping all the constant factors, gives 

TO > -y^, (22) 

as Qs 
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which is precisely the lower bound advocated by Arnold and Lenaghan in 

In general, using Eq. in Eq. to replace f{g^Nc)N^ with rj/T^ and 
dropping the numerical factors we conclude that 



This is the central result of this Section. It appears that the matching procedure 
between the CGC initial conditions and viscous hydrodynamics leads to depen- 
dence of thermalization time on shear viscosity. Moreover, the dependence shown 
in Eq. (|23l) implies that lower viscosity leads to shorter equilibration time! This 
conclusion agrees well with the fact that hydrodynamic simulations demand small 
shear viscosity and early thermalization time to describe RHIC data. It appears 
that, if our matching procedure captures correctly at least some features of the 
actual non-perturbative thermalization dynamics, short thermalization time and 
small viscosity of the quark-gluon plasma may be related. 

Finally, assuming that 't Hooft coupling g'^Nc is large at RHIC, one can use 
the result of [ 5 that f{g^ Nc) = tt/S in the strong coupling limit (indeed, of A/" = 4 
SYM theory), to estimate the thermalization time tq using Eq. H21() . Using as = 0.3, 
ce — Qs — 1-4 GeV for central RHIC collisions, and n{g) = 47.5 for three flavors 
yields tq ~ 0.5 fm, in a very good agreement with thermalization time required by 
hydrodynamic simulations at RHIC. 
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